Mem. S.A.It. Vol. 75, 282 
©SAIt 2004 



Memorie deiia "tj 




The Black Hole Mass of 
Abell 1836-BCG and Abell 3565-BCG 

o 

E. Dalla Bonta^'^, L. Ferrarese^, E. M. Corsini\ J. Miralda-Escude^ , L. Coccato^, 
55 . and A. Pizzella^ 

' Dipardmento di Astronomia, Universita degli Studi di Padova, Vicolo deirOsservatorio, 
I 3, 1-35122 Padova, Italy, e-mail: elena.dallabonta@unipd.it 

^ Herzberg Institute of Astrophysics, Victoria, Canada 
' Institut de Ciencies de I'Espai (CSIC-IEEC)/ICREA, Bellaterra, Spain 
(— I ' Max-Planck-Institut fuer extraterrestrische Physik, Garching bei Muenchen, Germany 

6 



C/3 



Abstract. Two brightest cluster galaxies (BCGs), namely Abell 1836-BCG and 
Abell 3565-BCG, were observed with the Advanced Camera for Surveys (ACS) and the 
Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope. 
By modeling the available photometric and kinematic data, it resulted that the mass of 
Abell 1836-BCG and Abell 3565-BCG are M. = 4.8!jJ^ x 10' Mq and M. = 1.31°;^ x 10' 
Mq at Icr confidence level, respectively. 
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1. Introduction lution of SBHs and the hierarchical build-up 

' of galaxies is most directly probed by the sys- 

O ; Many nearby galaxies have revealed large dark ^^^^ ^^^^ ^^^^ undergone the most extensive 

, masses in small regions, so far explained only protracted history of merging, namely the 

O ■ in terms of supermassive black holes (SBHs, ^^^^-^^ ^^^^^-^^ associated with the largest 

see for a review Ferrarese & Ford 2005). gg^s. This is the case of brightest cluster 

Spectroscopic and photometric data at high j^^^^^ ^^^^^y j^^^. j ^^^^^^ 1^. 

spatial resolution have inade it possible to de- i^inosities, as well as their having a merging 

H , rive relatively accurate SBH masses, M., for ^-^^^^^ ^^-^^ unmatched by galaxies in less 

a number of galaxies. These measurements ^^^^^^^ environments, make BCGs the most 

are rehable if the SBH sphere of influence is proj^ising hosts of the most massive SBHs in 

solved. It was also found that M. correlates ^j^^ j^^^j Universe. Therefore from Laine et al. 

with other galaxy properties, such as bulge (2003) it was selected a sample of BCGs, pre- 

luminosity (Magorrian et al. 1998, Marcom ^-^^^^ j^^^^ ^1^^ ^^^^ ^^^^-^^ ggj^^ ^^^^^^ 

& Hunt 2003), light concentration (Graham _ ^ _ relations), and 

et al. 2001), and bulge velocity dispersion cr ^^-^^ ^pj^^^^ influence could 

(Ferrarese & Merritt 2000; Gebhardt et al. ^^^^j^^^ ^-^^ ^^-^ 

2000). The inextricable link between the evo- ^ Telescope Imaging Spectrograph (STIS) 



X 



Send ojfprint requests to: E. Dalla Bonta 



on board the Hubble Space Telescope (HST). 
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Fig. 1. Observed [N II] /!6583 kinematics (filled circles) along with the best-fitting model (solid line) for 
the SBH mass of Abell I836-BCG and Abell 3565-BCG. The observed and modeled velocity curve (top 
panels) and velocity dispersion radial profile (central panels) are shown for the slit along the major axis 
(central panels), and the two offset positions (left and right panels) of the gas disk. The corresponding 
modeled flux profile is shown in the bottom panels. 



2. Imaging and spectroscopy 

Abell 1836-BCG and Abell 3565-BCG were 
observed with three filters of the Advanced 
Camera for Surveys (ACS): F435W (which re- 
sembles the Johnson B filter), F625W (simi- 
lar to r in the Sloan Digital Sky Survey pho- 
tometric system), and the narrow -band (w 130 
A) ramp filter FR656N, covering the redshifted 
Ha and [N II] AA654S, 6583 emission lines in 
the nucleus of each galaxy. The images were 
used to determine the optical depth of the dust, 
stellar mass distribution near the nucleus, and 
intensity map. Besides, high-resolution spec- 
troscopy of the Ha and [N II] /16583 emission 



lines was obtained for both galaxies, to mea- 
sure the central ionized-gas kinematics. A ma- 
jor axis spectrum was obtained with the slit 
crossing the nucleus of the galaxy. Two addi- 
tional offset spectra were obtained displacing 
the slit by one slit width on either side, per- 
pendicularly to the slit axis. Ionized-gas kine- 
matics was measured by fitting position, width, 
and flux of the [N II] /16583 line. Figs. 1 and 
2 show the kinematics of Abell 1836-BCG 
and Abell 3565-BCG, respectively. They re- 
veal regular rotation curves and strong central 
velocity gradients, thus they are good candi- 
dates for dynamical modeling. 
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Fig. 2. distribution for Abell 1836-BCG (top 
panels) and Abell 3565-BCG {bottom panels) as 
a function of SBH mass {left panel), inclination 
{central panel), and stellar mass-to-light ratio {right 
panel). The dotted horizontal lines indicate the con- 
fidence levels on the best fitting values. 



3. Dynamical modeling 

For this study, the procedure described in 
Coccato et al. 2006, to which the reader is re- 
ferred for details, was followed to model the 
ionized-gas kinematics. Briefly, a synthetic ve- 
locity field was generated assuming that the 
ionized gas is moving in circular orbits in an 
infinitesimally thin disk centered at the nuclear 
location, under the combined gravitational po- 
tential of stars and SBH. The model is pro- 
jected onto the plane of the sky for a given 
inclination of the gaseous disk, and then de- 
graded to simulate the actual setup of the spec- 
troscopic observations. The latter step includes 
accounting for width and location (namely po- 
sition angle and offset with respect to the cen- 
ter) of each slit, STIS PSF, and charge bleed- 
ing between adjacent CCD pixels. The free pa- 
rameters of the model are the mass M, of the 
SBH, mass-to-Hght ratio {M/L)i, of the stellar 
component (and dark matter halo), and inclina- 
tion / of the gaseous disk. Both (M/L)* and / 
are assumed to be radially invariant. Although 
i can be estimated from the images of the disk, 
the ionized gas is concentrated in the inner- 
most region of the dust disk and slight warps 
can be present even in the case of disks with a 



regular elliptical outline (e.g., Ferrarese et al. 
1996; Sarzi et al. 2001; Shapiro et al. 2006). 
Therefore ; is best treated as a free parame- 
ter The surface brightness distribution of the 
ionized gas is treated as input, by the narrow- 
band imaging. The velocity dispersion of the 
emission lines was assumed to be described 
by the radial function cr{r) = ctq + crie^''^''', 
with the choice of the best parameters to re- 
produce the observables. M,, {M/L)i,, and / 
were determined by finding the model param- 
eters which produce the best match to the ob- 
served velocity curve, obtained by minimizing 
= Ziv-Vmodf/SHv) where v + 5(v) and v,„orf 
are the observed and the corresponding model 
velocity along the different slit positions, re- 
spectively. 

4. Results 

For Abell 1836-BCG a three-dimensional grid 
of models was explored, with < M, < 3.0 x 
10'° Mo in 2.0x10*^ Mo steps, 0° < ; < 90° 
in 1° steps, and < (M/L)* < 40 (M/L)o in 
0.5 (M/L)o steps. The best model fitting the ob- 
served rotation curve requires M, - 4.8^^^ x 
10^ Mo, / = 76 + 1°, and (M/L)* < 4.0 
(M/L)o, where the errors on M, and /, and 
the upper limit on {M/L)i,, are quoted at the 
Icr confidence level. This is the largest SBH 
mass to have been dynamically measured to- 
date. The model is compared to the observed 
[N II] /t6583 kinematics in Fig. [1] Figure |2] 
shows Icr, 2cr, and 3cr confidence levels indi- 
vidually on M,, i, and {M/L)^,, according to 
the Ax^ variations expected for one parame- 
ter, with the other two held fixed at their best- 
fitting values. 

The rotation curves of Abell 3565-BCG 
along the three slit positions were fitted for 
a grid of model parameters defined by < 
M. < 5.0 X lO** Mo in 2.0x10^ Mo steps, 
0° < < 90° in 1° steps, and < (M/L)* < 15 
(M/L)o in 0.4 (M/L)o steps. The best model re- 
quires M. = 1.3+J}^ x lO** Mo, / = 50 + 1° and 
(M/L)*= 9.0 + 0.8 (M/L)o, where all errors are 
given at the 1 cr confidence level. The model is 
compared to the observed kinematics in Fig.[T] 
Fig. |2] shows the confidence levels on M,, i, 
and {M/L)i, alone. 
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Fig. 3. Location of the SBHs masses of our BCGs 
(stars) with respect to the near-infrared M.-Lb,iige 
relation by Marconi & Hunt (2003). NGC 1399 was 
added (Houghton et al. 2006). The SBH masses 
based on resolved dynamical studies of ionized gas 
(open circles), water masers (open squares), and 
stars (filled circles) are plotted. Dotted lines repre- 
sent the Icr scatter of the relation. 



5. Conclusions 

Several authors have recently focused their at- 
tention on the position of BCGs and massive 
galaxies in the upper end of the SBH scaling 
relations, using both theoretical and, even if in- 
direct, observational considerations. Bernard! 
et al. (2007) found that BCGs define a shal- 
lower cr - X relation than the bulk of early- 
type galaxies. They interpreted it as due to 
a non-linear correlation between logA4. and 
logcr. This would bring the Al.-cr to underes- 
timate M, for high cr values. Furthermore, ac- 
cording to von der Linden et al. (2007) BCGs 
follow a steeper Faber-Jackson relation than 
non-BCGs. This implies that BCGs can follow 
at most one of the power-law relations, either 
the M,-cr or the M,-L. Both Bernard! et al. 
(2007) and von der Linden et al. (2007) used 
SDSS data. On the contrary, Batcheldor et al. 
(2007) argue that SBHs masses predicted from 
NIR luminosities are consistent with masses 



predicted from cr. They attributed the discrep- 
ancies to the presence of extended blue en- 
velopes around the BCGs. 

In order to compare the SBH mass determi- 
nations of Abell 1836-BCG and Abell 3565- 
BCG with the prediction of the near-infrared 
M,-Li,uige relation by Marconi & Hunt (2003) 
the 2MASS //-band luminosities of the two 
BCGs were retrieved from the NASA/IPAC 
Infrared Science Archive. Figure [3] shows the 
location of the SBH masses in the M,-Li,uige 
plane. Only the SBH mass of Abell 1836-BCG 
is not consistent with the relation, since it is 
larger than the one expected (i.e., 8.9x10^ M©). 
The same is true only for M 87, among the 
seven galaxies with M,> 10'' dynamically 
measured so far. The cr of Abell 1836-BCG is 
not available in literature, thus it is not possi- 
ble to discuss its behavior in the M.-cr rela- 
tion. On the other hand, all the galaxies with 
A[,> 10'' Mq lie on the relation. There is no 
strong evidence that the upper end of the scal- 
ing relations is regulated by a different law. 

The only way to hush-up the debate is 
adding new M,. Before the advent of new 
classes of telescopes from earth and space, 
adaptive optics is a viable solution (e.g., 
Houghton et al. 2006). It will be possible to 
understand the behavior of both At, -cr and 
M, -X relations and, in case, if BCGs and nor- 
mal galaxies have a bimodal trend. 
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